Correlation between Charge Inhomogeneities and Structure in Graphene and Other 

Electronic Crystalline Membranes 



Doron GazitQ 

Institute for Nuclear Theory, University of Washington, Box 351550, Seattle, WA 98195, USA 

(Dated: October 5, 2009) 

Only one atom thick and not inclined to lattice defects, graphene represents the ultimate crys- 
talline membrane. However, its structure reveals unique features not found in other crystalline 
membranes, in particular the existence of ripples with wavelength of 100 — 300 A. Here, I trace 
the origin of this difference to the free electrons in the membrane. The deformation energy of the 
lattice creates a coupling between charge fluctuations and the structure, resulting in ripples on the 
membrane, correlated with charge inhomogeneities. In graphene this mechanism reproduces the 
experimental result for both charge puddles and ripples. 

PACS numbers: 61.46.-w, 73.22.-f, 74.90.-(-n, 87.16.D-. 



Physical membranes are objects in which one of the di- 
mensions is small compared to the other two, giving them 
an effective two-dimensional (2D) character. An intrigu- 
ing class of physical membranes is that of the crystalline 
membranes, that have a solid structure, usually of a tri- 
angular or hexagonal 2D lattice [T]. 

Beautiful examples of such systems, that spread in 
magnitude and scale, exist in our world. In biology, a 
famous example is the cytoskeleton of red blood-cells, 
whose structure is vital for the operation and stability of 
the cell, that forms a triangular crystalline lattice built of 
spectrin proteins. In soft condensed matter, one can cre- 
ate crystalline lattices by polymerizing liquid interfaces 
[T]. Recently, an ingenious experimental method, using 
mechanical cleavage, has conquered the final limit, isolat- 
ing graphene - a single-layer of carbon atoms, organized 
in a hexagonal lattice [2]. The same method has since 
been used to isolate individual crystal planes of other 
layered materials |3]. 

With such a vast variety of crystalline membranes in 
nature, it is of no surprise that understanding the struc- 
tural properties of these systems has attracted the atten- 
tion of many physicists [ll|4l[5l|6l[a[8l|9l[10l[IIl[12]. The 
fact that graphene can be used to construct nanometer- 
sized electronic applications, has only enhanced the need 
of a profound understanding of its structure jl3j . This 
interest intensifies, in view of the Mermin- Wagner theo- 
rem, which forbids the existence of long range order in 2D 
systems due to diverging thermal vibrations. This seem- 
ing contradiction between experiment and theory, which 
peaked with the discovery that graphene is stable even 
when it is free standing, i.e., without the support of a 
substrate [2] , is resolved by allowing out-of-plane fluctu- 
ations, that induce frustration between the large thermal 
vibrations in 2D and the competing gain in elastic energy. 
This frustration stabilizes a globally flat phase at finite 
temperatures |T]. This "almost flat" phase is character- 
ized by a scale invariant structure, at wavelengths much 
longer than a characteristic size At = ,,3^" rj. , deter- 
mined by the bending energy k, and the 2D- Young mod- 



ulus Kq {T is the temperature) [TJHISllSlEKHllSllinilll]- 

However, experimental studies of graphene have revealed 
rather different features. 

Meyer et al [T^ have isolated a free-standing graphene, 
thus demonstrating its long-range-order and stability. 
Furthermore, they used transmission electron spec- 
troscopy to study its structure. Their finding, which 
was since reproduced by other experimental groups, is 
that the graphene sheet exhibits spontaneous rippling, 
with amplitude of about 3 — 10 A, and wavelength esti- 
mated to be A = 100 - 300 A [H [13 ES]. Lattice defects 
were not found, and thus cannot account for this non- 
vanishing curvature [HI [T7] . Considering the fact that 
in graphene At ~ 12 A, these ripples clearly violate the 
scale invariance that should govern at this scale. 

As graphene is the ultimate crystalline membrane, this 
difference has to be addressed theoretically. In this Rapid 
Communication, I suggest that the ripples in graphene 
are a signature of the fact that it is not a regular crys- 
talline membrane, since it has an additional degree of 
freedom - the free electrons that occupy its tt band ~ 
thus it is a representative of a different class of materials 
- electronic crystalline membranes, in which an interplay 
exists between the electronic and structural degrees of 
freedom [T8j [19] . This interplay leads to the excitation 
of ripples in electronic crystalline membranes in general, 
and graphene in particular. 

To reach this conclusion, I start by modeling the 
"almost-flat" phase of a membrane. In-plane deforma- 
tions are characterized by a two dimensional vector fleld 
u, and out-of-plane deformations by a fleld h. When 
considering the equilibrium state of the electronic crys- 
talline membrane, without allowing charge fluctuations 
in the conduction electrons, the mesoscopic structure of 
the membrane can be described by thermal fluctuations 
around this equilibrium using the elastic free energy: 

F[u, h]^^-J d'x [^{V'hf + 2^iul + Xul] , (1) 

where Uij is the strain tensor, Uij = ^ {diUj + djUi) + 



2 



\{dih){djh) (summation over repeated indices is implied 
throughout the manuscript). The coefficients are the 2D 
elastic properties of the membrane. In graphene, ex- 
periments have verified this approximation, as no lat- 
tice defects were found even at strains ~ 10% [131 [T7] . 
and the elastic constants were estimated: a bending en- 
ergy K « l.leV, bulk modulus A -|- /Lt ~ 7.3eV A , and 
shear modulus /i « 5.7 eV A pO] (the resulting Young 

modulus is Kq — "^a^^+A^^ ~ 13eVA ). /it and A were 
estimated from the sound velocities. These elastic prop- 
erties originate in the a band, which is a consequence of 
the in-plane sp^ hybridization, that forms a deep valence 
band, and partly by the tt band, which is perpendicular 
to the plane. This was also verified by Monte-Carlo sim- 
ulations with realistic, though phenomenological, inter- 
atomic potentials (that cannot take into account charge 
fluctuations) [2T1[2HI23]. An additional veriflcation that 
Eq. ([ij indeed describes such simulations, is due to the 
fact that the structure of graphene predicted by them ex- 
actly fits the theory of crystalline membranes, quantita- 
tively predicting the scale invariance and the anomalous 
exponents of the bending energy and elastic constants at 
long wavelengths [23], however unable to reproduce the 
ripples. Thus, the ripples have a different origin. 

The free electrons, which in graphene occupy the half- 
filled TT band, differentiate an electronic crystalline mem- 
brane from regular crystalline membranes. Indeed, the 
free electrons couple to the structure through a defor- 
mation energy. The source of this deformation potential 
is the local change in the Fermi energy measured from 
the bottom of the valence band, proportional to the local 
change in area 5S. The resulting deformation potential 
has the form Vs — = Duu {a is the nearest neighbor 
distance), where D is the Fermi energy of the 2D elec- 
tron gas [24'. In graphene this deformation energy was 
found to be the main source of deviations of the elec- 
trical transport properties from ballistic transport. The 
specific value of D is a matter of debate. However, using 
charge carrier mobility measurements of electron doped 
graphene, one achieves a value oi D k, 29 eV, which com- 
pares well with other estimates of this energy [25 . 

In the presence of spatially varying density of 
TT-electrons 5n{x), the deformation energy is just 
J cPxDuii{x)dn(x). [32] . 

The structure of the electronic membrane is thus de- 
termined by the following free energy: 



F[u, h, Sti] ^ SeJH + ^ J (fxK{\/^h)^ 



(2) 



+ - J (fx [2fiufj + Xu% + 2DuuSn\ . 

Where Eee is the electron-electron energy due to the 
charge density. The free energy evidently couples be- 
tween between elastic deformations and charge inhomo- 
geneities. 



In order to understand the behavior of this system, let 
us first concentrate in the electron-electron interaction. 
Characterizing this interaction is not a trivial task, as it 
includes a solution of a strongly correlated many body 
problem. However, one can estimate it by: 
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with e the charge of the electron, the static dielectric 
constant £ (originates from screening of the electron- 
electron interaction), and Sn{q) a Fourier transform of 
6n{x). The value of the static dielectric constant is not 
fully known, since the effective fine structure constant in 
graphene is of order unity. Recent perturbative analysis 
of the electron-electron screening by Kotov et al. [5B] 
has shown that the static dielectric constant is £ « 3 — 4. 
However, they demonstrated that the perturbative anal- 
ysis receives large corrections at higher orders of per- 
turbation theory, concluding that additional screening is 
expected in the non-perturbative solution. Here, I will 
take this value as a starting point for the analysis, and 
examine the effect of increasing the value of the static 
dielectric constant on the final conclusions. 

The theory is now quadratic in the charge fiuctuations, 
thus can be written as: 



F[u, h, Sn] 
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where A(g) = A — §^1- The theory is gaussian also in 
the in-plane deformation field u, thus these degrees of 
freedom can be integrated out, resulting in an effective 
free energy that depends only on the out-of-plane defor- 
mations: 



Feff[h] = - 



d'^q 
(2^ 



(5) 



Where <^[h{x)] = ^(5, 



effective Young modulus, given by: 



^^)dihdjh and K(q) is an 
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Clearly, for long wavelengths (corresponding to g — > 
0) K{q) Kq, implying that long wavelength sound 
waves can be used to measure the elastic constants, as 
they appear in Eq. ([T]) , neglecting the effects of possible 
formation of charge inhomogeneities. 

In order to investigate the possibility of ripple excita- 
tions, one is interested in the behavior of the electronic 
crystalline membrane at finite wavelengths. Observing 
Eq. ([6]), it is clear that interesting phenomena occur 
around the length scale = = e'^(2ti+\) • -'^^ particu- 
lar, for qo > q > the effective Young modulus is 
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negative, thus representing a true competition with the 
bending energy K{V^hy, allowing height fluctuations. In 
order to quantify this, one has to solve the thermodynam- 
ics of the system described by Eq. (|5| . 

In general, Dyson equations can be written for the 
scale evolution of the effective bending rigidity Kn^q) = 
and Young modulus KR{q): 



GCq.TjATo 



1 /qt 

2 q 



(7) 
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Here, qx 



KokBT 



E((7) is the sum of all IPI two-point 



diagrams, and ^'(q) is the sum of all IPI four-point di- 
agrams. These equations determine the structure of the 
electronic crystalline membrane. Evidently, two intrinsic 
length scales exist, i.e., qo, originating in the electronic 
degrees of freedom, and qr, in which thermal effects be- 
come significant. The existence of the scale qo differenti- 
ates electronic crystalline membranes from regular crys- 
talline membranes, whose structure is controlled only by 
the thermal scale. In the current case, however, there is 
an interplay between the two length scales. 

As aforementioned, K{q) Kq for q <^ q^. Hence, 
in this limit, the Dyson equations reduce to the usual 
case of regular crystalline membranes, extensively stud- 
ied particularly in the thermodynamic limit L ~f oo 
(g ^ 0) (see, e.g., [Illl|5J[i|71[Hl[a[Tni[IIJ [T5]). As 
a result, electronic crystalline membranes, for which this 
hmit corresponds to the regime q <C {qo, Qt}, inherit the 
same behavior at long wavelengths, i.e., have a stable, 
asymptotically flat, phase, with the IPI diagrams gov- 
erned by anomalous exponents, viz. ^ q^"'' and 
q^"''", with 77 « 0.8 and = 2 - 2ry. 

However, the behavior outside this regime, and espe- 
cially about the electronic scale qo, is governed by the 
special functional structure of the Young modulus, cf. 
Eq. (|6|. In particular, Eq. ([s]) shows that the formal 
structure of the effective Young modulus survives the 
scale evolution, i.e., there always exists a region of nega- 
tivity for the young modulus, with a zero at g = 20^1o 
and a singularity around q ~ qo- As a result, the physics 
at g ~ (7o will be governed by this length scale, rather 
than qT- Due to this, as well as the fact that finite q 
behavior is of interest, we calculate S(g) in the one-loop 
approximation (V denotes principal value): 



S(g) = V 



(Pk K{qk) \q X 



(9) 
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the normal-normal corre- 
^ kBT/{KR{q)q^) = 
as it indicates enhanced cor- 



We search for maximum 
lation function: G{q) 

{K/Ko)[{q/qT? + n<l)]\ 
relation, which will manifest itself as ripples on the mem- 
brane [21 , if the maximum is pronounced, i.e., if its width 
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FIG. 1: (Color online) The normal-normal correlation func- 
tion as a function of the wave-number q for different temper- 
atures. 
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FIG. 2: (Color online) Ripples wavelength as a function of 

^(2m+A) ■ 



the parameter = ^ ^ 



is smaller than the characteristic wavelength. The am- 
plitude of the ripples is determined by the maximal value 
of the correlation. 

The resulting dimensionless normal-normal correlation 
function KoG{q) / n, as shown in Fig. [l]for various tem- 
peratures, presents a rather pronounced peak. The main 
temperature dependence is in qT, since simulations sug- 
gest rather weak temperature dependence of the elas- 
tic constants and the bending rigidity [52]. Evidently, 
the temperature dependence of G{q) vanishes for q « 
0.02 A ^ . This signatures the onset of the long-wavelegth 
regime, where the one-loop approximation given here is 
expected to fail. Comparing the value of the correlation 
function at this wavelength to its maximal value is a mea- 
sure of the amplitude of the ripples. It can be seen that 
the ripples are more pronounced for higher temperatures, 
and have shorter wavelength. The width of the correla- 
tion maximum weakly depends on the temperature, as 
it is mainly affected by the electronic length scale qo, 
rather than the thermal one. The temperature depen- 
dence should be further investigated, taking into account 
changes in other effects, e.g., defect formation and elastic 
constants. 
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holds the largest uncertainty in the model. 



due to the static dielectric constant. Thus, in Fig. [2] 
the wavelength maximizing the correlation function is 
plotted as a function of (at room temperature) , for a 
physical range of values of this parameter, chosen around 
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^o{D = 29eV,£ = 4) « 20 A. As seen in the plot, the 
maximal correlation occurs for wavelengths in the range 
^ = 100 - 300 A. This wavelength region reproduces 
the ripples found in experiments on suspended graphene 

Thus, one expects a stable, asymptotically flat, phase 
which exhibits charge fluctuations and ripples, with the 
same correlation length. The character of this correlation 
can be analyzed through Eqs. (4][5), where the charge 
fluctuations and in-plane deformations were integrated 
out. The average values of these integrated out observ- 
ables are readily recovered, and a relation between the av- 
erage charge fluctuation and the average lattice deforma- 
tion is established (Sn) cx {S[h{x)]), where S[h] = — V^<i> 
is the Gaussian curvature of the surface. This is a key 
difference between the current work and previous theo- 
retical studies of the correlation between charge puddles 
and height fluctuations, e.g., Ref. [18], that found that 
charge fluctuations are proportional to the mean curva- 
ture, i.e., to V^/i (indeed, a recent experiment ^TT did 
not find such a correlation) . 

In suspended graphene, no experiment has studied 
both the topology and the charge inhomogeneity to- 
gether. However, the correlation lengths of the two disor- 
der phenomena were measured in different experiments, 
finding ripples with characteristic size of 100 - 300 A 
[m [T5j [16], and a correlation length of 300 A for charge 
puddles [2H]. In the latter experiment, though graphene 
was suspended on top of Si02, it was shown that the 
substrate had no effect on the structure of the charge 
puddles. Experiments that probed both phenomena si- 
multaneously were accomplished only for graphene on 
top of Si02 substrate [ST] [29], but have shown signifi- 
cant substrate effects. It is important to note that in a 
setting where the graphene is located very close to the 
substrate, its structure and its charge fluctuations would 
be pinned to the surface structure and impurities on the 
substrate, and not as discussed here. In addition, the cur- 
rent work discusses neutral graphene. Doping is analogue 
to external stress whose sign corresponds to the majority 
charge carriers, thus decreasing (increasing) the rippling 
in the case of electron (hole) doping [H]. This effect 
can explain the results of Ref. [29] , where the graphene 
sheet was doped by an external gate voltage. A different 
external source for ripples in graphene is adsorption of 
molecules [TSJ [30] • The current work is different since it 
proposes the electrons as an intrinsic source for ripples 
at thermodynamic equilibrium. 

In conclusion, charge puddles and ripples in graphene 
are found to be a signature of the fact that graphene is 
not a regular crystalline membrane, but the herald of a 
new class of materials - electronic crystalline membranes, 
demonstrating strong interplay between the dynamics of 
the free electrons in the membrane and its mesoscopic 
structure. Clearly, this implies that had the tt elec- 
trons not been free, the ripples would vanish. Indeed, 



graphane, an insulating graphene derivative in which 
each carbon atom is connected to a hydrogen atom, was 
found to exhibit reduced corrugations [3T] . 

This paper offers a theoretical approach to charac- 
terize the two main intrinsic disorder phenomena in 
graphene, i.e., charge inhomogeneity and structural de- 
formations. As graphene is promising material for tech- 
nological use, understanding disorder phenomena and 
correlations among them is essential for a successful de- 
sign and quality control of future applications. 

I thank K. Novoselov, L. Radzihovsky, Dam T. Son, 
A. Fasolino and G. Bertsch for helpful discussions. This 
work was supported, in part, by DOE under grant no. 
DE-FG02-00ER41132. 



* Electronic address: |doron.gazit@mail.huji.acH| 
[1] D. Nelson, T. Piran, and S. Weinberg, eds.. Statistical 
Mechanics of Membranes and Surfaces, Second Edition 
(World Scientific Pubfishing, 2004). 
K. S. Novoselov et al., Science 306, 666 (2004). 
K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. 
Khotkevich, S. V. Morozov, and A. K. Geim, Proc. Natl. 
Acad. Sci. 102, 10451 (2005), arXiv:cond-mat/Q503533. 
M. J. Bowick and A. Travesset, Phys. Rep. 344, 255 
(2001). 

Nelson, D. R, and Peliti, L., J. de Phys. 48, 1085 (1987). 
F. David and E. Guitter, Europhys. Lett. 5, 709 (1988). 
Guitter, E. , David, F. , Leibler, S. , and Peliti, L. , J. de 
Phys. 50, 1787 (1989). 

J. A. Aronovitz and T. C. Lubensky, Phys. Rev. Lett. 
60, 2634 (1988). 

P. Le Doussal and L. Radzihovsky, Phys. Rev. Lett. 69, 
1209 (1992). 

C. Schmidt, K. Svoboda, N. Lei, L Petsche, L. Berman, 

C. Safinya, and G. Grest, Science 259, 952 (1993), 
http: / /www. sciencemag.org/cgi/reprint/259 / 5097/952. pdf. 
Mark J. Bowick, Simon M. Catterall, Marco Falcioni, 
Gudmar Thorleifsson, and Konstantinos N. Anagnos- 
topoulos, J. de Phys. I 6, 1321 (1996). 

D. Gazit, Phys. Rev. E [in press], arXiv: 0907.3718. 
A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

J. C. Meyer et al., Nature 446, 60 (2007). 
V. W. Brar, Y. Zhang, Y. Yayon, T. Ohta, J. L. McCh- 
esney, A. Bostwick, E. Rotenberg, K. Horn, and M. F. 
Crommie, Appl. Phys. Lett. 91, 122102 (2007). 
V. Geringer, M. Liebmann, T. Echtermeyer, S. Runte, 
M. Schmidt, R. Riickamp, M. C. Lemme, and M. Mor- 
genstern, Phys. Rev. Lett. 102, 076102 (2009). 
Y. Zhang et al., Nature 438, 201 (2005). 

E. -A. Kim and A. H. Castro Neto, Europhys. Lett. 84, 
57007 (2008), arXiv:cond-mat/0702562. 
D. Gazit, Phys. Rev. B 79, 113411 (pages 4) (2009), URL 
http : //link . aps . org/abstract/PRB/v79/ell3411 

F. Guinea, B. Horovitz, and P. L. Doussal, 
Solid State Communications 149, 1140 (2009), 
ISSN 0038-1098, recent Progress in Graphene 
Studies, URL |http: //www. sciencedirectT] 



[9 
[10 

[11 



[12; 

[13 

[14 
[15 



[16 

[17 
[18 

[19 

[20 



I c om/scieiice/article/B6TVW-4VWBlBJ- 6/2/ | 
I ff 31a08bc6090791b54393d5bca0d72e 
[21] A. Fasolino, J. H. Los, and M. I. Katsnelson, Nature 

Materials 6, 858 (2007). 
[22] K. V. Zakharchenko, M. I. Katsnelson, and A. Fasolino, 

Phys. Rev. Lett. 102, 046808 (2009). 
[23] J. H. Los, M. L Katsnelson, O. V. Yazyev, K. V. Za- 
kharchenko, and A. Fasolino, Phys. Rev. B 80, 121405 
(pages 4) (2009), URL jlittp : //link, aps . org/abstract/, 
I P RB/v80/e 1 214051 

[24] H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 
(2002). 

[25] K. L Bolotin et al., Phys. Rev. Lett. 101, 096802 (2008). 

[26] V. N. Kotov, B. Uchoa, and A. H. Castro Neto, Phys. 
Rev. B 78, 035119 (2008). 

[27] A. Deshpande, W. Bao, F. Miao, C. N. Lau, and B. J. 
LeRoy, Phys. Rev. B 79, 205411 (pages 6) (2009), URL 
http : //link . aps . org/abstract/PRB/v79/e205411 

[28] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. 
Smet, K. von Klitzing, and A. Yacoby, Nature Physics 
4, 144 (2008), 0705.2180. 

[29] Y. Zhang, V. W. Brar, C. Girit, A. Zettl, and M. F. 
Crommie, ArXiv e-prints (2009), 0902.4793. 

[30] R. C. Thompson-Flagg, M. J. B. Moura, and M. Marder, 
EPL (Europhysics Letters) 85, 46002 (4pp) (2009), URL 
http : //stacks . iop . org/0295-5075/85/46002 

[31] DTX;. Elias, R. R. Nair, f. M. GTMohiuddin, S. V. 
Morozov, P. Blake, M. P. Halsall, A. C. Ferrari, D. W. 
Boukhvalov, M. I. Katsnelson, A. K. Geim, et al.. Sci- 
ence 323, 610 (2009), (in the supplementary material), 
http: / /www. sciencemag.org/cgi/reprint /323 / 5914/610. pdf, 
URL http : //www, sciencemag. org/cgi/content/l 

I abstract/ 323/5914/610 

[32] Note that charge neutrality demands J (fx5n{x) = 0. 



